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Abstract 

A Particle Tracking Velocimetry experiment has been performed in a turbulent flow at interme- 
diate Reynolds number. We present experimentally obtained stretching rates for particle pairs in 
the inertial range. When compensated by a characteristic time scale for coarse-grained strain we 
observe constant stretching. This indicates that the process of material line stretching taking place 
in the viscous subrange has its counterpart in the inertial subrange. We investigate both forwards 
and backwards dispersion. We find a faster backwards stretching and relate it to the problem of 
relative dispersion and its time asymmetry. 



I. INTRODUCTION 



The time asymmetry existing in a turbulent flow governed by the Navier-Stokes equation 
has recently been discussed in the connection with the diffusion properties of two nearby 
fluid particles initially close to each other. A stochastic model study by Sawford et al. |l8l | 
showed that backwards dispersion is faster than the corresponding forwards dispersion. This 
has practical implications for a variety of applications dealing with transport and mixing in 
turbulent flows. 

The asymmetry was experimentally verified by a Particle Tracking Velocimetry (PTV) 

n 

experiment by Berg et al. [lj who speculated that kinematic infinitesimal material line 
stretching has its counter-part in the inertial range. Different stretching rates could then 
explain the observed difference in forwards and backwards dispersion. Through a finite 
Reynolds number formulation of the much celebrated Richardson law (r 2 ) = get 3 it was 
found that following two particles backwards in time the dispersion was faster by a factor of 
approximately two compared with following the particles forwards in time. This ratio was 
also obtained from Direct IN umerical Simulation (DNS) analysis by the same authors (data 
originally presented in 

HQ) 

. Sawford et al. [18| found the same ratio to be much larger: 

between 5 and 20. 

Yet another stochastic model was developed by Luthi et al. Q|. By assuming that 
stretching rates behave self-similar results consistent with [1] were obtained. By self-similar 
stretching we think of particle separation which is similar on the different scales. The self- 
similarity of stretching rates has, however, not been shown. 

In this paper we will explore forwards/backwards dispersion in the context of coarse- 
grained strain dynamics. The latter has recently got a lot of attention in the turbulence 
commu nity since it is more or less related to Large Eddy Simulations (LES) 



community 
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One can expect Kolmogorov similarity scaling to hold for the coarse-grained 
quantities as long as the filtering scale is less than the integral scale in the flow where the 
turbulence is solely represented by the kinetic energy dissipation of the flow, e 
As in we define the coarse-grained velocity gradient tensor to be 

du A 
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where u A (x, t) is the coarse-grained velocity over scale A: 

5f(x) = i/ dru,(x + r). (2) 

Here -Ba/2 denotes a Ball with radius A/2 and its volume. 

We will follow the lines of Q and Luthi et al. and present experimental evidence of 
self-similar stretching in a turbulent flow of intermediate Reynolds number. The result is 
linked to particle pair separation and is able to account for the difference in forwards and 
backwards dispersion. 



II. EXPERIMENTAL METHOD 



We have performed a Particle Tracking Velocimetry (PTV) experiment in an intermediate 
Reynolds number turbulent flow. The flow has earlier been reported in PTVJjs an 



experimental method suitable for obtaining Lagrangian statistics in turbulent flows 
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16l l22j : Lagrangian trajectories of fluid particles in water are obtained by tracking 
neutrally buoyant particles in space and time. The flow is generated by eight rotating 
propellers, which change their rotational direction in fixed intervals in order to suppress a 
mean flow, placed in the corners of a tank with dimensions 32 x 32 x 50cm 3 (see FigQ). The 
data acquisition system consists of four commercial CCD cameras with a maximum frame 
rate of 50Hz at 1000 x 1000 pixels. The measuring volume covers roughly (12cm) 3 . We use 
polystyrene particles with size ~ 400/zm and density very close to one. We follow (9(1000) 
particles at each time step with a position accuracy of 0.05 pixels corresponding to less than 
lOyum. 

The Stokes number, tj/t^ (tj denotes the inertial relaxation time for the particle to the 
flow while t v is the Kolmogorov time) is much less than one and the particles can therefore 
be treated as passive tracers in the flow. The particles are illuminated by a 250W flash 
lamp. 

The mathematical algorithms for translating two dimensional image coordinates from the 
four camera chips into a full set of three dimensional trajectories in time involve several cru- 
cial steps: fitting gaussian profiles to the 2d images, stereo matching (line of sight crossings) 
with a two media (water-air) optical model and construction of 3d t raj ectories in time by 
using the kinematic principle of minimum change in acceleration 15J, |2l( . 
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FIG. 1: Experimental setup 
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TABLE I: Turbulence characteristics obtained from the von Karman model 14]: Mean kinetic energy 
dissipation e, Kolmogorov length r\ = (i^ 3 /e) 1 / 4 , Kolmogorov time t v = (v/e) 1 / 2 , integral length Li nt , 
integral time Tl , velocity fluctuations o\ = | (cr^ +&u y + a u z ) ■ The Reynolds number is defined as Re\ = ^f- 
with the Taylor micro scale A = ' '"' rT ~ 



From eqn. ^we can define the coarse grained strain (s 2 ) and (Q 2 ) where (sy) and (ttij) 
are the symmetric and anti-symmetric part of the coarse-grained velocity gradient tensor 
(Aij). The eigenvalues and eigenvector of (s^} is denoted (Aj) and (Ai) respectively. Due to 
incompressibility ^j(Aj) = 0. In addition we define (Ai) > (A 2 ) > (A3) so that the most 
positive principal direction of strain is (Ax). 

The filtered velocity field (eqn. |2J) is approximated by least square fits of spherical in- 
compressible and orthogonal linear polynomials to discrete velocities of fluid particles inside 
spherical balls with diameter A: 



mm 



fl(A/2) 



d x AijXj^j 



(3) 



At least four particles are necessary to describe a three-dimensional shape and hence to 
estimate (Ay). Four particles form a tetrahedron and is the backbone of the analysis by 
Chertkov et al. [7]. In [10] we show that obtained from only four particles is quite far 
from the definition in eqn. |2]and therefore that using only four particles are not sufficient. 

We find that at least twelve particles are needed in order to obtain a reasonable ap- 
proximation of eqn. El of the coarse-grained quantities. Using this many particles has the 
drawback that we can not study the dynamics at the smallest scales. In Figure 121 (a) the 
radial distribution of particles is shown. The probability Np(r) of finding a given number 
iV of particles on the surface of a ball centered in our measuring volume with radius r is 
observed to increase with r 2 up to ~ 200^. This means that the particle density can be 
considered uniform up to this scale. For lager radius the density drops down because of non- 
perfect illumination at the boundary of the measuring volume. The cumulative distribution 
is show in (b) and is interesting because it gives us an estimate of the number of particles 
we can expect to find in a ball with radius r. The number is, however, only an upper bound 
since the ball is centered in the measuring volume and not around all individual particles in 
the flow. 




FIG. 2: Radial distribution of particles where N is the total number of particles in the volume, a) The 
probability of finding a particle in a spherical shell with thickness dr and with distance r from the center of 
the measuring volume, b) Cumulative distribution of (a). 

The smallest scale for which (u) and hence (s 2 ) and (Q 2 ) can be resolved is 8O77. A higher 
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particle seeding density which again makes particle tracking more difficult can decrease this 
number. All results reported use a minimum of twelve particles. 

As already reported in ^, the mean flow is slightly straining though with a charac- 
teristic time scale many times larger than the integral time scale of the flow. An alternative 
way of emphasizing the large scale mean strain is shown in Figure EH cos 2 (xj, Ai), where Xi 
is the three coordinate axes, is observed to peak for large r in the vertical direction while 
the horizontal axes decrease in agreement with an axis-symmetric flow. A slow convergence 
towards isotropy is observed as r = A is reduced. 
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FIG. 3: Squared cosine between the most stretching principle coarse grained eigenvector, Ai, and the three 
coordinate directions. Blue is the vertical axis-symmetric direction while the green and red is the horizontal 
directions. 



III. COARSE-GRAINED STATISTICS 

In the viscous subrange the stretching rate of infinitesimal material lines governs particle 
pair separation. The stretching rate is defined as 
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where d/dt denotes Lagrangian differentiation (following the particles). In the viscous sub- 
range Lit) when rescaled with the Kolmogorov time scale r„ becomes constant after a short 
time . In the viscous sub range the second order Eulerian longitudinal structure 

function f{r) = ({[Mj(x + r) — Uj(x)]r;/r} 2 ) is given by 



f(r) 



<f!>, 
15 ' 



r <^ rj 



(5) 



Through the definitions e = 2u(s 2 ) and r 2 = vje we can form r 2 = r 7(30/(r)). Thus 
motivated by viscous subrange scaling we move to the inertial range and define a time scale 



**(r) = 

In the limit r <C ri we have t+ = V2t„. 



2r 2 



15/(r) 



r < L 



(6) 
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FIG. 4: Coarse grained strain (s 2 ) (long dashed blue) and (f2 2 ) (long dashed red) as a function of coarse 
graining scale r. The solid black line is r -4 / 3 . The dotted lines are (s 2 ) and (f2 2 ) multiplied by i 2 . The 
vertical black dotted line indicates the integral length scale. 



As a function of scale r = A we plot (s 2 ) and (fl 2 ) in FigureEJ Both quantities are in the 
inertial range observed to be in agreement with the Kolmogorov similarity prediction r _4//3 . 
For r -> the ratio (fi 2 )/(§ 2 ) -> 2. 
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Also in Figure |U we have plotted (s 2 ) and (Q 2 ) multiplied with t 2 . In the inertial range 
we find that (s 2 )^ ~ 1 and (fi 2 )t 2 = 1.85. It thus seems that as defined in eqn. serves 
as a characteristic time scale for (s 2 ) and therefore t* ~ 1/ a/ (s 2 ). 

In [n| _/"(?") in eqn. IHlwas defined through a model valid for homogeneous and isotropic 
turbulence. In this paper f(r) are obtained directly from data. The behavior therefore is 
therefore slightly different. 



IV. PARTICLE PAIR SEPARATION 

The rescaled eigenvalues (Aj)t* are shown in Figure |SJ The trademark of turbulence, 
namely (A2) > which is necessary for both positive mean enstrophy and strain production 
is observed for all values of r. A slight decrease in (A2) is observed as r is increased which 
could be taken as a sign of the coarse-graining field becoming more Gaussian and hence 
(A 2 ) = 0. Both (Ai)t* and (A 3 )t* are constant and so are their ratio |(A3)|/|(Ai)|. For a 
direct comparison with viscous result we have divided the results with 
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FIG. 5: Eigenvalues (Aj)i* of (sij). i = 1: green, i = 2: blue and i — 3: red. The black curve is the ratio 
I (A3) l/l (Ai) I- 




When time is running forward two particles in a mean strain field will on average separate 
from each other along the direction of (Ai). In the backward case they will separate along 
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the direction of (A3). Since (A 2 ) > and therefore |(A 3 )| > |(Ai)| backwards separation is 
the faster one. This is shown in Figure El for times ts < t < 5ts where ts = { r o/^) 1 ^ is the 
Batchelor time, characterizing the time for which the initial separation should be regarded 
an important parameter in the separation process p . By closer inspection of Figure H we 
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FIG. 6: Angular dependence between (Aj) and the separation vector r for times ts ^> t 3> 5ts- Red curves 
are forwards dispersion while blue are backwards dispersion. 

can see that in the backward case r is slightly more aligned with (A3) than it is aligned 
with (Ai) in the forward case. This is due to the positiveness of (A 2 ) and a corresponding 
increase in the alignment of r with (A 2 ) in the forward case compared to in the backward 
case is therefore also observed. The values of cos 2 (r, Aj) are in the forward case 0.50, 0.32 
and 0.19 for i = 1,2,3 respectively. In the backward case the same values are 0.19, 0.28 and 
0.53. 

Stretching rates rescaled by i* are presented in Figured as a function of t/t B for different 
initial separation of pairs. 

For times t ~ 0.5tg the forward stretching and backward stretching occurs with the same 
speed. The separation vector is still randomly oriented and has therefore not yet aligned itself 
with the principal directions of the strain field. For longer times the backwards stretching 
rates saturates at 0.235 ± 0.005 while the forwards saturate at 0.18 ± 0.005. Whereas all 



9 



0.4 



0.3 




1 2 3 4 5 

t/t B 

FIG. 7: Stretching rates L(t) rescaled by t*. The blue curves correspond to forwards separation while the 
red curves corresponds to backwards separation (time running backwards) for different initial separations 
tq. Particle pairs with ro from 3 — 4mm to 9 — 10mm in steps of 1mm are shown. We have divided the 
y-axis with V2 in order to match the viscous limit. 



curves in the forward case collapse for times up to 3 — the curves in the backward 
case do not collapse so nicely. Because backward separation is faster the particles leave the 
measurement volume earlier: the larger r , the earlier. The particle pairs we thus observe 
for large times are likely to be slow pairs with low stretching rate. To reduce the effect of 
finite volume we only choose pairs which start within a small sub volume (30mm) of the full 
measuring volume. There is, however, no systematic way in which we can totally neglect 
the effect of a finite volume. 

The maximum stretching rates would occur if the particle separation vector r is fully 
aligned with the principal strain axes. In this case the rescaled stretching rates would 
be (Ai)i* and —(A^t* in Figure El for the forward and backward case respectively. The 
forward value Lf{t)/\/2/t± = 0.18 ± 0.005 is close to the values obtained in the viscous sub 
range flSS 3. 

Berg et al. |l( hypothesized that, if the particle separation was perfectly aligned with the 
principal strain axes, the ratio between forwards and backwards dispersion rates character- 
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ized by the Richardson- Obukhov constant ratio gb/gf could be determined as |(A 3 )|/ (Ax)) 3 . 
From Figure lathis number equals 1.13 3 = 1.44. From the calculated stretching rates we find 
that this is certainly not the case since 0.235/0.18 = 1.31. With the same argumentation 
this would give a ratio gb/gf = 2.23 ± 0.33 within errors of what was measured directly in 
Q. Although L{t) f < | (Ax) | and L(t) b < |(A 3 )| we thus have L(t) b /L(t) f < |(A 3 )/(A 1 )| and 
hence a larger ratio than in the fully aligned case. It is important to stress that the argu- 
mentation of how to relate the ratio gb/gf to L(t)b/ L(t) f through the Richardson-Obukhov 
law is not a rigourously derived result but merely a hand waving argument. 

V. CONCLUSIONS 

We have given evidence of the existence of, what we call self-similar stretching rates: 
when scaled by a relevant time t± scale which is a function of the second order structure 
function stretching rates of infinitesimal material lines has its counterpart in the inertial 
range. Furthermore it turns out that this relevant time scale t± is related to the coarse- 
grained strain field. The stretching is like in the viscous range far from being perfectly 
aligned with the most positive direction of strain which we have shown would lower the 
ratio between forwards and backwards dispersion. 

Whether or not the Lagrangian stretching rates found in this paper are universal or 
specific for this particular turbulent flow other experiments and/or DNS will have to decide. 
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